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Abstract 

Within the framework of the scalar-tensor theory (STT), its second post-Newtonian (2PN) ap- 
proximation is obtained with Chandrasekhar's approach. By focusing on an iV-point-masses system 
as the first step, we reduce the metric to its 2PN form for light propagation. Unlike previous works, 
at 2PN order, we abandon the hierarchized hypothesis and do not assume two parametrized post- 
Newtonian (PPN) parameters 7 and f3 to be unity. We find that although there exist 7 and f3 
in the 2PN metric, only 7 appears in the 2PN equations of light. As a simple example for appli- 

m 

cations, a gauge-invariant angle between the directions of two incoming photons for a differential 
f — , ■ measurement is investigated after the light trajectory is solved in a static and spherically sym- 

o: 

CN ■ metric spacetime. It shows the deviation from the general relativity (GR) <5#stt does not depend 

on /3 even at 2PN level in this circumstance, which is consistent with previous results. A more 



^ ■ complicated application is light deflection in a 2-point-masses system. We consider a case that 

H ! 

the light propagation time is much less than the time scale of its orbital motion and thus treat 
it as a static system. The 2-body effect at 2PN level originating from relaxing the hierarchized 
hypothesis is calculated. Our analysis shows the 2PN 2-body effect in the Solar System is one 
order of magnitude less than future ~ 1 nas experiments, while this effect could be comparable 
with 1PN component of <5#stt in a binary system with two Sun-like stars and separation by ~ 0.1 
AU if an experiment would be able to measure 7 — 1 down to ~ 10~ 6 . 
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I. INTRODUCTION 



Some future space missions, such as the Laser Astrometric Test Of Relativity (LATOR) 
1, 2\, the Phobos Laser Ranging (PLR) [3j, the Beyond Einstein Advanced Coherent Optical 
Network (BEACON) ^1, the Telemetrie InterPlanetaire Optique (TIPO) fl, the Astrody- 
namical Space Test of Relativity using Optical Devices (A.STROD) |6| and the Search for 
Anomalous Gravitation using Atomic Sensors (SAGAS) |7|, will measure distances of laser 
links and angles among these links with unprecedented precision. As a sensitive and useful 
tool in gravitational physics, especially for some high order effects, the propagation of light 
carries lots of information about the nature of spacetime and plays an important role in 
high-precision experiments and measurements. 

Thus, from the practical and theoretical aspects, it motivates us to investigate a 2PN 
light propagation model at the c -4 level in the framework of the scalar-tensor theory for a 
gravitational A-point-masses system. 



A. The necessity of 2PN c 4 terms 
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For the precision of distance measurements, TIPO could achieve millimeter level |5l and 

n n 

BEACON could be even higher, reaching 0.1 nanometer(nm) [4(. As explained in Ref. [8] 
(in Sec. V), a complete metric of the Solar System up to c -4 level is demanded for modeling 
the light propagation in those experiments. For the angle measurements, the precision of 
LATOR could achieve nano-arcsecond (nas) level or higher l|, |2j. It would require c -4 terms 
in the metric to calculate the deflection of light ray, because when light grazes solar limb 
the 2PN deflection contributed by the Sun is G 2 m^ ) /{c 4L R^ ) ) ~ 10~ 12 ~ micro- arcsecond 
(/zas), which is much larger than the threshold of LATOR. The corrections of alternative 
theories of gravity to this effect at 2PN order would be at least several orders of magnitude 
less than ~ /ias. However, these corrections perhaps still need to be considered for future 
experiments. 

Several authors have obtained the 2PN metric for the general relativity (GR) and alter- 
native theories of gravity. Chandrasekhar and Nutku first calculated the 2PN metric and 
equations of hydrodynamics in GR js]. A field theory approach is also employed to de- 
rive the 2PN Lagrangian for the equations of motion of iV bodies in the multi-scalar-tensor 
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theory without solving the metric lOj. With introducing an intermediate-range gravity 
term, the 2PN approximation of the scalar-tensor theory was obtained [ll] in which the 
energy-momentum tensor was expressed by using the invariant density [l2[. The 2PN ap- 
proximation of Einstein-aether theory [13| was deduced in the form of both superpotentials 
and an iV-point-masses. More recently, the IAU2000 resolutions are extended to include all 
the c -4 terms for the requirements from some space missions [8] . Most of such works express 
the 2PN metric in terms of superpotentials without definition of the masses and multipole 
moments of the bodies which is not trivial in the post-Newtonian order especially for the 
self-gravitating bodies in an iV-massive-bodies system. 

An explicit application of them is to model light propagation in the near zone of a grav- 
itational system into the level of the next leading order. The 2PN light deflection for one 



jody in the standard parametrized post-Newtonian (PPN) formalism 14j, [l5| was studied 



la, LLZ( by introducing another parameter (e or A) at c 4 of g^j, the space-space component 



of metric, under the isotropic gauge. The 2PN light propagation in the Schwarzschild space- 
time was detailedly researched in three gaug es (standard, harmonic and isotropic) through 
introducing two coordinate parameters 18]. Ref. [3] discussed the significance of some 
experiments to observe 2PN light effects for one body in the isotropic gauge. A practical 
relativistic model for the 2PN light propagation was developed in the harmonic gauge of 



general relativity 20], in which the 1PN contributions from the rectilinear and uniform N- 
body and the 2PN contributions only from the Sun were considered. The light propagation 
in 2PN framework of a stationary gravitational field for the Schwarzschild metric (one body) 
was formulated in harmonic gauge with introducing one parameter e at spatial isotropic and 
anisotropic terms of c for g^ [2JJ. The 2PN deflection of light in a spherically symmetric 
body was discussed in GR in three gauges (standard, harmonic and isotropic) [22j. All of 
above works mentioned are only considering one body in the 2PN order. The 2PN gravita- 



tional redshift in GR was derived 23] by superpotentials. Under the f(R) theory, the 2PN 



weak lensing was explored in the isotropic gauge 



24|. 



B. The reason for a scalar-tensor theory 

Although Einstein's general relativity has passed nearly all the tests in the Solar System, 
alternative theories of gravity are still required for deeper understanding of the nature of 
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spacetime and for testing possible violations of the Einstein equivalence principle (EEP) in 



the forthcoming more precision level 



251 ] . In order to testing and distinguishing alternative 



gravitational theories, the PPN formalism introduces ten parameters in a post-Newtonian 



metric to include various gravity theories 



15(. However, the PPN formalism is only 



restricted to the 1PN approximation. Some authors (e.g. 26H28j ) discussed how to param- 
eterize the 2PN metric. To extend the PPN formalism, one possible way might be to derive 
2PN metric for various gravity theories and then find the independent superpotentials in 
the 2PN level. After that, parameters in the 2PN order could be introduced and endowed 
with some meanings. However, in this paper, we only focus on one gravitational theory and 
discuss its parameters. 

Among these alternative theories, the most eminent case is a scalar-tensor theory (STT) 
because it is the simplest and most natural way to modify GR. Many modern theories, such 
as the extra-dimensional theory, the string theory, the braneworld and the noncommutative 
geometry, which try to unify gravity and microscopic physics or to explain the dark energy in 
cosmology, demand a scalar field besides the metric tensor (see 29[ for a review). Especially, 
Damour and Nordtvedt proposed [30] that the deviation of the PPN parameter 7 from 1, 
presenting the discrepancy between STT and GR, might range from ~ 10" 7 to ~ 10" 5 . 
The contribution of this deviation in light deflection at 2PN order could be as large as 
(7 — 1)G 2 7TIq/(c 4 Rq) ~ 10~ 16 ~ 0.02 nas, which might be important for future detection as 
well. Therefore, we shall work in this paper with the scalar-tensor theory of gravity. 

Motivated by forthcoming space experiments involving propagation of light in the Solar 
System, some researchers had studied 2PN light propagation in STT. In the case of a one- 
body system, 2PN light deflection and light propagation were reported in Refs. 3jj and 



32J . As a more comprehensive work, the scalar-tensor propagation of light was investigated 



in Ref. 



33| , in which it neglected the differences of PPN parameters 7 and (5 from 1 at c 



level according to present experiment results and used the hierarchized hypothesis so that 
only Sun's contribution remained at c -4 level while influences from Sun and planets were all 
included at c~ 2 and c~ 3 terms. 

To make our results fully based on STT, we will keep 7 and /3 at c -4 level which requires 
a PPN definition of mass, such as Eq. (!33|) . 
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C. The extension to an iV-point-masses system 



As metioned in Ref. 



33], the hierarchized hypothesis has its own limitation. In the 



practical point of view, by abandoning this hypothesis which means extending to an iV-body 
system, the solution of light equations might improve the accuracy which could be achieved 
33[ | . As a rude estimation, one coupling term in 2PN metric Yl a Sfe^ a G 2 m a m\ ) l [c^r a Tb\ 



where r a = \x — y a \, = \x — and the trajectories of the a-th and 6-th masses are 
respectively represented by y a (t) and yb(t), could reach ~ 10~ 16 (equivalent to ~ 0.1 nas), 
if we consider the Solar System (a hierarchic system) that a-th mass is the Sun and the 6-th 
mass is Jupiter when the light grazes Jupiter's limb. If the influence to the background light 
sources by a binary system with comparable masses are considered, the coupling term in 
2PN might raise larger contributions. Extension to an iV-body system would also be helpful 
to determine which terms have to be included for a specific mission and could be a good 



33]. 



test-bed to evaluate the accuracy of the hierarchized model 

In the theoretical point of view, it is a natural development to build a 2PN theory of a 
gravitational iV-body system for its dynamics and light propagation within it which might 
show some subtle but interesting effects due to the non-linearity in the 2PN order. In 
principle, the 2PN light propagation under the STT should include relativistic multipolar 
moments of each body in the iV-body system. 1PN global and local metrics with the 
definitions of the multipolar moments and spins in STT have been given in Refs. HH 



Therefore, in this paper, we consider 2PN light propagation in an iV-point-masses system 
under STT as our first step. 



In what follows, our conventions and notations generally follow those of Ref. [36]. The 
metric signature is (—,+,+, +). Greek indices take the values from to 3, while Latin 
indices take the values from 1 to 3. A comma denotes a partial derivative, and dot over 
a quantity denotes a derivative with respect to time. Bold letters denote spatial vectors. 
The plot of this paper is as follows. In the Sec. II, the 2PN metric and equation of light 
in the scalar-tensor theory are given. Subsequently, in Sec. Ill, we reduce the results of 
Sec. II to a 2PN metric of an A-point-masses system for light propagation by following the 
method used in Ref. 37|. And the 2PN light equation is obtained in this condition. In 



Sec. IV, applications will be given. In Sec. IIV A\ we derive the 2PN light-ray trajectory 
and deflection in a static, spherically symmetric spacetime and study the parameters in 2PN 
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terms by comparison among the scalar-tensor theory and others. In Sec. IIVP4 we calculate 
the light deflection angle caused by a 2-point-masses system and estimate the magnitudes 
of effects for different cases. Finally, the conclusion and discussion are outlined in Sec. V. 



II. 2PN METRIC IN THE SCALAR-TENSOR THEORY AND EQUATIONS OF 
LIGHT 



A. 2PN metric in the scalar-tensor theory 



The action for the scalar-tensor theory based on 34| reads 
,3 r / Q(X\ i 6n 



s 



16tt 



(1) 



where c is the speed of light, 9(<p) is an arbitrary coupling function of the scalar field </>. 
R and g = det(g M „) denote respectively the Ricci scalar and the determinant of the metric 
tensor g^ v . The matter field is denoted by From Eq. (GQ), we can see the matter fields 
\l/ only interact with the metric field (namely, g^v). This means the trajectory of a free-fall 
test particle only depends on the spacetime geometry so that it satisfies EEP. Although 
violations of EEP at galactic and cosmological scales can not be ruled out, we still focus on 
the scalar-tensor theory satisfied EEP in this paper. 
Variation of the action (JT]) with respect to g a ^ has 



R 



, 2 (t^ 2 g,„T 



(2) 



where D fi 
by 



)\iiv9 iiv ' i T^ u is the stress-energy-momentum tensor of matter which is defined 
_ d(^Cj) d d(y/=gd) 



c 
~2 



-gT f 



[IV 



dg 



dx c 



dg 



flU 



(3) 



and T is the trace of T^ v . Following Refs. [38|, |39[, the mass, current, and stress densities 
can be defined as 



cT°\ 



a 



an = c 2 TK 



(4) 
(5) 
(6) 
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Variation of the action ([IJ with respect to yields 



Based on Ref. 



34|, the scalar field is expanded as its background value 0o as follows 

= 0o(l + C), (8) 

where ( is a dimensionless perturbation of the background value O - Especially, decomposi- 
tion of the coupling function 9((j)) reads 

9{<t>) =u + wiC + ^ 2 C 2 + • • • , (9) 

where Co> = #(0o) ; w i = (d^/dC)^^,, • • ., uj n = (d n 9 /d^ 71 )^^, which can lead to the gauge- 
invariant definitions of PPN parameters 7 and (3 (see Eqs. ( 1A3I) and ( 1A10I) in Appendix 
A) 

u + l 



7 = 



(3 = 1 + 



w + 2' 



(2w + 3)(2w + 4) 2 ' 

By using Chandrasekhar's approach [j], Q], we deal with the theory in the form of a 
Taylor expansion in the parameter e = 1/c. The expansions of the metric gr^ and the scalar 
field to the second order have the forms of 

000 = -l + e 2 N + e*L + e 6 Q + 0(7), (10) 
g 0i = e 3 Li + e 5 Qi + 0(6), (11) 
9lj = b i3 + e 2 E ij + e A Q^ + 0(5), (12) 

(2) (4) (6) 

C = e 2 ( +e 4 C +e 6 ( + (9(7), (13) 

where 0(n) means of the order e n . 

According to Eqs. flS}, © and (ITS]) , we derive 

(2) /I (2) 2 (4)\ 

Wo + e 2 WlC +W_ W2C +WlC j + o(5) ; (14) 
^lL + ^ + 0(4)). (15) 
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We use the gauge condition imposed on the component of the metric tensor proposed by 
Kopeikin & Vlasov as follows: 




Based on fields equations of Eqs. <$fy and ©, we obtain the evolution equations of the 
metric coefficients at 2PN order by using the gauge conditions Eq. ffT6l) (see Appendix A 
for detail). 



B. 2PN equations of light 



Generally, for a photon propagating in a spacetime in which Einstein Equivalence Prin- 
ciple (EEP) is valid, the basic equations of light based on Ref. 4l| are 

dx^ dx u 



= g^- 



d 2 x l 



e r°, 



dt dt ' 

dx % 

va 



t \ dx v dx a 

va J 



dt 2 \" V ° dt ' V °) dt dt 

where we have replaced affine parameter with coordinate time t. Following Ref. 
suming x = csfi and fj. ■ fi — 1, we find the expression for s from Eq. (fl71) 

s = 1 - \e 2 (N + Hijti'fjP) - e 3 L kf i k 

+ l -e A [-Q l3 ^^ -L+ l -NH l3 ^^ - l -N 2 + ^H^^H^/j, 1 }. 



(17) 
(18) 



2l|, as- 



(19) 



Then, by substituting the metric Eqs. (TTUT) - (TT2T) into Eq. ( Tl8|) . we obtain the equations of 
light propagation based on Eq. ( fl8l) as follows 



x? x k 



• H • k 

X X 



X J X 



\ 1 -' -' r 

o^V + n^jk,i ^ij,k A „• 

2 2 c c c c c c 



+£ 



Li j h L 



X 



■3 



X 



- „JV, t - + -H ik , 



3,1 c 2 x,,t c 2" iM c c c 



H, 



x J 



ij,t ' 



-L 



Sjq3 X^ 

c c c 



X 3 X k 1 



X^ X 



+£ 2 { — -^H ik N k — NNa h H u Hij k -H u Hj k l 

'2 cc c c 2 cc 



x 1 x 3 



x 3 x k 1 



x^ x^ 



'Li t + A — L j Qij.k 1" nQjk.i f + 

2 cc cc 2 c c 1 



(20) 



From Eq. ( 120]) . we can see that the 2PN metric for light propagation could be cut off to 



£00 = -1 +e 2 N + e 4 L + 0(5), 



(21) 



(22) 
(23) 



g 0i = e 3 Li + 0(5), 
g tj = 5ij + e 2 E ij + e^Q i3 + 0(5). 

Furthermore, with using the relationship: Hij = 5ij^N based on Eqs. (1A5I) and (1A6I) in 
Appendix A and by substituting x ■ x = c 2 s 2 and Eq. f[T$j) into Eq. f[2"Uj) . the equation of 
light is simplified as 
1 



-(l + 7 )^-(H- 7 )JV, fc -- 
2 c c 



+ L 



x 



3,i o V 1 + !f)N,t n L j,k o-^fcj 

c 2 c 2 c c c 2 c c e 



1 Qlj^ 1 •37"^ tJC^" 



+^<j _ 7 (i + ^nn. _ L . t + V + ( 7 2 _ _ ^Q ijtk — — 



l 



V x k 1 



T.Qik,j 1- 7:Qjk.i L j > + 0(3), 

2 cc 2 c c c c 1 



III. 2PN LIGHT-RAY PROPAGATION IN A SYSTEM OF N POINT MASSES 
A. Solving the metric for 2PN light propagation 

Considering an A^-body system of nonspinning point masses as our first move, we follow 
the notation adopted by Ref. |37| and use the matter stress energy tensor as follows: 

c 2 T^(x,t) = J2»a(t)vZv v Ax - y a (t)), (24) 

a 

where S denotes the three-dimensional Dirac distribution, the trajectory of the a-th mass 
is represented by y a (t), the coordinate velocity of the a-th body is v a = dy a (t)/dt and 
v% = (c, v a ), and \i a denotes an effective time-dependent mass of the a-th body defined by 



v a v a 



(25) 



where subscript a denotes evaluation at the a-th body and m a is the constant Schwarzschild 
mass. Another useful notation is 

,2" 



1 + ^ 

c z 



(26) 



where v 2 = v a ■ v a . Both \i a and \i a reduce to the Schwarzschild mass at Newtonian order: 
Ha = rn a + 0(2) and jl a = m a + 0(2). Then the mass, current, and stress densities in Eqs. 
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JSJ) and (jSJ) for the N point masses read 

o" = ^2^ a 5(x - y a {t)), 

a 
a 



(27) 
(28) 
(29) 



(2) 



The next step is to work out N, C an d in 1PN approximation as 
N = 2n- 1 {-AnGa} 

= 2 Y. G {y + ^(/i a r a )} + 0(3) 



Gm a ^ ( Gm a 



+EE 



G 2 m a m b 



Av 2 a - (n a v a f 



n a n ab ) } + 0(3), 



+ 2 ( 2 -37)EE 



a 



G 2 m a m b 
r a r a b 



a b^a 



ab 



(30) 



( ^ (1 _ 7)E ^ +£2 /I (1 _ 7)E 



4i^ - (n a t> a )' 



+ o-#-CE ^ + 5(1 - •/) E E ^(-w) 



ab 



+0(3), 



(31) 



K - (n a v a f 



+ 2(2"37)EE 



G 2 m a mb 



a b^a 



r a r a b 



EE 

a by^a 



G 2 m a m b 



ab 



n a n ab ) \ + 0(3), 



(32) 



by the relation of /x a that 



/x = m a < 1+e 



m a { 1 + e' 



(2-37)X;^ + k' 

17^ ^ab 2 



+ 0(4) 
+ C(4)h 



(33) 



where r a = \x — y a \ and r ab = \y a — Vb\- Scalar products are denoted with parentheses such 
as (n a v a ) = n a ■ v a , n a = (x- y a )/r a and n ab = (y a - y b )/r ab . For L u 

U = -2(1 + 7 )D-'{-4^} = W G,n f f _> 7 A = £ + 0(2).(34) 

J \x z\ v a 
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Based on Eq. (1A8I) . L can be simplified as 

L — (37 — 2/3 — 1)D- 1 {-4ttG(7} + 2(7 - l)D- 1 {-47rG'a fefc } - i/3iV 2 



/3N 



(3 7 - 2/3 - 1) G{^(N) a \ + 2(7 - 1) £ g{^} - 5 

'a 'a , , 1 a 1 ab 

a a a b^a 

, , ' a' 

a 



with the help of 



'a a' ab 

a a b^a 



And 



i/3iV 2 = - 2 ^^-2^Vy^ + 0(2). 
2 V rl r a r b 



C = 2(7 -1)> u a + (1 + 67 -37 -4/3 > > 

'a , , ' 0' ab 

a a b^pa 

a a a b^a a 



-0(1). 



It follows that we derive the metric Qij as follows 



DQij = -8(1 + 7 )7r6?(7 - - ^(1 + 7)^^ 



+^|87rG<7(^7 - \i 2 - /3 + l)iV + 16 7 7rGa fefc + ^( 7 2 +/3 - l)V 2 (iV 2 ) 
where the quadratic part of potentials in Qj,- can be rewritten as 



= 5 E G!m «( S S + % A > (4) + 4 E E CWAA; (^) 



a 



12 



where d a i denotes the partial derivative with respect to y a . The integral of the self-terms 
can be readily deduced from Alnr a = l/r„ {3?]]. On the other hand, the interaction terms 
are obtained by 



AlnS, 



1 



ab 



r a r b 



(42) 



where S ab = r a + r b + r ab [12| . Consequently we solve 
-i(l+ 7 )D- 1 [iV iJ iV J ] 



- \(l + 7) Y, G2m - { d l ^ a 4)- 2 (!+7)EE G 2 m a m b d ai d bj In S, 



^(1 + 7)E G 



2_2 / n a n fe _ $ij 



C2 



2 ( 1 + 7)^^G ,2 m a m b 



r a bS, 



ab 



(43) 



where two relations that 



d% In r a 



dij - 2rt a n{ 



(44) 



and 



OaiObj in ^ a fe = h 



q2 



(45) 



are used. Then, can be worked out as 



Qii = 2(1 + 7) E ^« + 5(1 + 7) E ^« - 2(1 + 7) E E 



G 2 m a m b 

r abS a b 



-2(1 +7) EE 



G 2 m a m b 



a b^a 



s 



< b -oKb+ n i) 



al> 



+5i 



47 £ - 2(3 7 - 3 7 2 - 2/3 + 2) £ £ G ' 2/ "" / " / ' 



+ (27 2 -i7 + 2/3-|)E^ + 2(l + 7 )5:$: 



, 'aTab 
a b^a 

G 2 m a m b 



a bj^a 



r abS, 



ab 



+2( 7 2 + /3-l)££ 



a b^a 



G 2 m a m b 



+0(1). 



(46) 



13 



B. The 2PN metric of light-ray propagation for an iV-point-masses case 



Collecting all these results together, we have the 2PN metric for light propagation in the 
scalar-tensor theory as 



fi'oo 



vl- 



-i + ^£^ + W- 2 /?£^ + 2 (i + 7)£^ 

a v a " a 

+2( 2 - 2f>- 1) y. E ^ + E E -2^EE G2m ° mt 



a b^a 



Vab 



a b^a 



a b^a 



»V6 



+0(5), 

-e 3 2(l + 7 )£^< + 0(5), 



(47) 
(48) 



^ + ^7 E ^ + { - ^7 + 2/3 - §) £ ^ - 7 E 

a v a " a 

+ E E i 2 ^ +/?-!)— + 2(2/3 - 7 - 2)— + 7 - (rW) 



a b^a 



Vab 



+ 2(1 + 7) 



r abS, 



ab 



+- 4 U (1 + 7) E ^<< + 2d + 7) E ^« 



-2(1 + 7 ) E G2 ™^ 

a b^a 



r abS, 



ab 



C2 
°ab 



+ 0(5). 



(49) 



If the above result returns to GR ( 7 = /3 = 1) and two-body case, it coincides with the result 



of Ref. 



37j. It is worth noting that the parameters at the spatial isotropic (^ XL ° ™ a ) 



|) and 



and anisotropic (J? G 7° r a r a) P ar ts of £ 4 f° r ^« are respectively (2 7 2 — ^7 + 2/3 

'a 

|(1 + 7) in the scalar-tensor theory (STT). For one-body case, different theories impose 
different values on the above parameters (see Table [J). When we take Einstein's general 

n n 

relativity (GR), the parameters coming from two parts are totally unity [20]. Ref. [13| 
shows, for a single body, Einstein-aether theory (AE-theory) has two unequal coefficients of 
the spatial isotropic and anisotropic part in the 2PN but only one parameter C14 presents 
in them. Recently, Ref. [21] introduces only one parameter e to parametrize these two parts. 
From our calculations, it shows the parametrization with single parameter e in Ref. 2jJ is 
not valid for the STT. Parametrized 2PN issues are more complicated and were researched 
to a certain extent by the previous works [26I+28I] . but it is beyond the territory of this paper. 
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TABLE I: Parameters at the spatial isotropic and anisotropic parts of e for in different theories. 



Parameter 


GR 


AE-theory 


Parameter in 


STT 


at 


[20] 


[13] 


[21] 


Our work 


r G 2 m 2 
°ij r ' 2 


1 


1 + |ci4 


e 


2 7 2 -i 7 + 2/3-| 


G 2 m 2 i j 
r 4 


1 


1 - 


e 


3(1 + 7) 



C. The light-ray propagation equations 



Substituting metric coefficients N for Eq. (1301) . for Eq. ( 152)) . Lj for Eq. (j34j) . L for 
Eq. ()35D and Q^- for Eq. flU into Eq. flU]), we have 



— -F^p^y + F^ PN + i 7 ^, 



(50) 



where the 1PN monopole component is 



1PN 



"(1+7)E 



n! - 2- 



(51) 



the 2PN monopole component is 



n a x 



2PN 



e2(1+7) £^J 2(1+7)<+ ^ 



.2 a 



ra„ — 



G 2 m a m b f 1 



+e 2 (i + 7)E E _^ n * _ _ KnabK + (i _ 7) _^ 

^ ^ r- I r a r a6 2r z 



a6 



1 



+ (l + 3 7 )-<-^< 6 - 



2 (n a x)x J (n a n a6 ) (n a x)x* 1 {n ab x)x l 



+ 



af) 



ab 



r lb c 2 



+2( i_ 7 )^(!^_ 2 (i_ 7 ). 1 



r b cr 



r a r b c z 



G 2 m a m b \ 1 (n afe x)(n a x) + (n afe x) (n fe x) 



a 6 



-n 



Tab 



ab 



-4- 



1 [(n ab x) + (n b x)][(n a x) + (n b x)] 



S, 



ab 



< b - K) + 



1 (n ab x) + (n b x) 



1 

-2— 



n\ b - <) + 2— 

Tab 



[n a x)x % {n b x)x l 



1 (™afcA)' 



(n a x) i x l 
— ifi 1 

c a c 

K + <\ 



, 1 [Kix) - (w a g)][(w a6 a;) + (n 6 x)] . 4 1 (n ab x) - (w a i;) 

+ ^~ ^ l n a + n bJ H 



06 



5, 
' ab 



(n b x) , x l 
(52) 
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the influence of the bodies' orbital motions is 
K = _ £2 (1 + 7) £ ~ —< + 



1 (v a x)(n a x) 
^{n a v a ) + - 2 



+^ 2 (l + 7)£ 



Gm a 



2 3, \ 2 i (^a^-j 



+ 



2t- 



2\ lb a 



2 l^a^ 



V n .X 



+ 2(n a v a ) — 
c c 



(53) 



In the case of a one-body system studied by many previous works, we always can construct a 
reference frame within which the body is static so that the velocity-dependent terms vanish. 

IV. APPLICATIONS IN AN iV-POINT-MASSES SYSTEM: N = 1 AND N = 2 



After obtaining the light-ray propagation equations, we will apply them into a 1-point- 
mass system (N = 1) and a 2-point-masses system (N = 2) to calculate the deflection angles 
for differential measurements. 



A. Deflection angle by a 1-point-mass system 



For a 1-point-mass system, it has a static and spherically symmetric spacetime. The 2PN 
light deflection angle caused by the central body has been computed in some previous works, 
such as Ref. 16l-l22l|. Here, we leave the details of calculation in Appendix [B] and show the 
final result of the deflection angle A# = 9 — -$o as 

1 + 7\ 4Gm a 



A9 



2 ) c 2 d 
(l + 7) f 



2(l + 7 )--(7 + 8 7 )vr 



G 2 ml 



(54) 



where d a = k x (r a x k) is an impact parameter to represent the closest distance between 
the unperturbed light ray and body a, d a = \d a \. When 7 = 1, Eq. (154")) will reduce to the 



result of GR 



4l|. And it is identical with the deflection angle reported in Ref. [32], which 



considers 2PN light propagation in a one-body system under the framework of STT. 

Table HT1 lists the summary of 2PN parameters appearing in the metric, light-ray trajectory 
and deflection angle in (1) GR; (2) AE-theory; (3) Ref.[2l|; (4) STT (our work). In GR, there 
is no any parameter. For AE-theory, there is only one parameter C14. From Ref. 21(, three 



16 



TABLE II: Summary of 2PN parameters in different theories. 



Theories 


Parameters 

ILL UlCUXll^ 


Parameters in 

li orrrf"— ts^it - f 1 f»r''l"ri'r'\?' 


Parameters in 

9PN rlpflprtinn 


GR[22] 


None 


None 


-8 + fvr 


AE-theory [13] 


C14 


C14 


-8+ + §c u V 


From Ref.[21] 




7,/3,e 


-2(l + 7 ) 2 + [2(l + 7)-/3 + f e ]7r 


STT [our work] 


7,/3 


7 


-2(l + 7) 2 + i(l + 7)(7 + 8 7 )7r 



parameters, 7 , (3 and e, show up. It is very interesting that only one parameter 7 appears in 
2PN light-ray trajectory and light deflection for the scalar-tensor theory although there are 
two parameters 7 and (3 in 2PN metric for the theory. If we let 7 = 7 — 1, Eq. (151]) returns 
to 

AO = 56 GR + 56 STT , (55) 



where 



r/i 4Gm a 15 G 2 m 2 

m a \ f R Q 



2 



1.75 — — arcsecond 

M Q )\d a 



2 / n \ 2 



+3 - 5 °I^Jllfl (56) 



_2Gm a _ /n 31 N G 2 m 2 , ^,_ 2 



2.1 x lCr 5 J \M Q J V d, 

+0(f). (57) 

Here, <5#gr is the light deflection equivalently caused by GR and its 1PN and 2PN effects are 
~ 1.75 as and ~ 3.50 /xas respectively, both of which are much greater than the thresholds 
of future ~ 1 nas space missions, 5#stt represents the discrepancy between STT and GR 
in the deflection angle. For the effect of the deviation from GR in STT in the deflection 
angle which would be measured by experiments conducted in the solar system, the leading 
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effect is 18.3 //as at the 1PN level. As Damour and Nordtvedt mentioned in Ref. 30], the 
deviation of 7 from unity at the levels of ~ 10~ 7 to ~ 10 -5 might be detected by some new 
space projects (e.g. LATOR). The 2PN term in 5#stt might also be important because it 
could reach ~ 0.1 nas, close to the thresholds of future space-borne experiments. 



B. Deflection angle by a 2-point-masses system 

From a 1-point-mass system to an iV-point-masses system, the spacetime will no longer be 
spherical symmetric and static. This complicates the problem significantly. We suppose that 
there are three characteristic length scales: the length scale of the iV-point-masses system £, 
the distance between the light source and the system D\\ and the distance from the observer 
and the system $K b s ; and two characteristic time scales: the time scale of orbital motions of 
the iV-point-masses system X and the time scale of the light crossing the system T cross . In 
general, £, 9i\ and SH b s could be arbitrary, while X is determined by the total mass of the 
system VJl and £, as well as X cr0 ss < £/c. 

In what follows, we will consider a specific case that the gravitational system contains 
only two point masses with body a and b (N = 2) and £ ^> m&x(Dl\, £H bs), and then 
calculate the resulting deflection angle. 



1. Light equations 

For this 2-point-masses system with the configuration £ 3> max^A, 5^obs), it indicates 
Across X. However, this does not mean that the impact of the orbital motion can be 
neglected, if this effect firstly emerges at less than or equal to 2PN level. Nevertheless, the 
calculation including this motion demands a huge amount of work and, thus, we would leave 
it for our next steps. In this paper, we firstly consider a case of this binary system that the 
effects of velocity-dependent terms are smaller than those of 2PN static effects, i.e. Ff, = 
in Eq. (1501 . so that the light equations are 

x l = F{ PN + FI pn \ibe + F^pn^bE) (58) 
where the contribution at 1PN is 

(n a x)x l 



F{ PN = -(1 + 7)— r- 



< - 2 , 



+ a <r+ b, (59) 
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the contribution of one-body effect at 2PN is 



G 2 m 2 



2PN\1BE 



n„x 



e\l + 7)— ^ 2(1 + T K + ^-nj, - 



{n a x)x 



(60) 



and the contribution of two-body effect at 2PN is 

Fl PN \2BE 



G 2 m a m b 



2r 



^{n a n ab )n l a + (1 - 7)r~r< 



ab 



r a r b 



+ (l + 3 7 )^- 



-n„ h - 



o r 2 "'ab 2 
a b I a' oh L 



2 (n a x)x l (n a n ab ) (n a x)x l 1 (n ab x)x % 

' O o o o 



oft 



ab 



r a r b c z 



+^ 2 (l + 7) 



G 2 m a mb 



1 (n ah x)(n (1 i;) + (n ab x)(n b x) 



-n 



Tab 



-4 



1 [(n a6 x) + (n b x)} [(n a x) + (n&i;)] 



S, 



n„ 



ab 



-2- 



,n b x 



1 



n 



a6 "a 



ni)+2 



1 



. ab ~ O + 

n a x)x l 



ab 



1 (n afe x) + (n b x) 



r a 
n b x)x l 



[n„,x 



+ 



1 {n ab xf 

Tab C 2 



X 

—K 

c c 



K + n ] 



, Q 1 [(ngtg) - (nqg)][(wafcg) + (n b a;)] . , ^ 1 (n a6 :r) - (n a x) 

n \ n a + n b) 



s, 



ab 



n 



[n b x ) , or 

~n\ - - 

c c 



-a-B-fc, (61) 
where the symbol a 6 means the same terms but with the labels a and b exchanged. 



K + ™ 



,2. Light-ray trajectory 

Following the procedure mentioned in Appendix [Bj we could work out the time derivative 
of light-ray trajectory -8x(t) and compute the deflection angle caused by these two point 
masses. For F^ PN and i^pjvli-BEj we could obtained corresponding trajectory of the light-ray 
by adopting the iterative method used in Ref. [la ]. However, in order to integrate FI pn \2bei 
we assume the light source, light-ray trajectory and the observer are all closer to the body 
a than to the body b, i.e. r a < r ab along the light path from end to end. Thus, we could 
obtain 

-x(t) = k + -5x 1PN (x N ) + -5x 2PN {x N )\ 1BE + -5x 2 p N (x N )\ 2 BE, (62) 

c c c c 

where 

1 . . . , Gm a [ k x (r a x k) ,1 . , Gm b [ k x (r b x k) \ N 

-5x 1PN (x) = - 1 + 7 )- r ±\ 4 } - + k - 1 + 7 )- T ±< 4 ~ + k ^ 3 

c c J r a [ r a - k ■ r a J c 2 r b { r b - k ■ r b } 
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-Sx 2 pn{x)\ibe 

c 



1 m ^G 2 m 2 a/1 , 1., ^G 2 m 2 
4( 1 +7)^f(fc-^K- i (l + 7)^f 

1 + 7) 2 



m 



G 2 m 2 
+ — r^d, 
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— + 



r a (r a - k-r a ) \r a r a - k ■ r a 



+ 



>+7)(7 + 8 7 )l 
G 2 m 2 



k- 



a 1 / 7T k ■ r a 

+_/ +arctan . 



4 



c 4 r a 



^fc 



i( 1 + 7)(5 + 4 7 )l-(l + 7 ) 2 -—l— 

1 'a ' a "» ' o 



-I(l + 7 )(7 + 87)^ 
cH b 



' ! (1+ 1 ) 
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1 / 7r k ■ r b 

+ 7 - + arctan — - — 
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){l + 7 )(5 + 47)- - (1 + 1? \— 

4 r b n-k- r b 



(64) 
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-{k ■ r ab )ln — 
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+ (l + 7) 2 
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ah 
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— + arctan • 

Tab V 2 



c 4 r ab 



k ■ r. 



l(l + 4 7 )ln ra + f- r " 



(65) 



where r 0a = r a {t = t ) and r 0a — |r n J. For GR (7 = = 1) and one body case, Eq. fl62|) 
will reduce to the results of Ref. [18]. 



5. Lig/ii deflection 

With the definition of a gauge-invariant angle # between the directions of two incoming 
photons (light signals 1 and 2) and we use that the position of the photon at the moment t 
of observation coincides with the position of the observer so that 

x bs = x 01 + c(t - to^kx + 8x x = x 02 + c(t - t 02 )k 2 + Sx 2 , (66) 

where (t i,*oi) denotes the moment and position of the light signal 1 of emission and 
(t 2, ^02) for the light signal 2 respectively. 

By assuming that the initial unit vector of light-ray 1 k\ is perpendicular to the line 
connecting the body a and the b, r ab , and light-ray 2 moves along the line connecting the 
body a and the observer, which lead to that geometrical relationships that 

ki ■ r a k 2 -r a fei • r ab d la 

« cos^o, = 1, = 0, — « sm(l (67) 

r a r a r ab r a 

d 2a = 0, 2 , la = sin^o, 2 Vab = sin$ , k 2 = — . (68) 
dia r ab r a 

and letting d\ a = d a , we could obtain the deflection angle for a static observer as 
A6 = 6 - ^0 

. /, >Gm a / cos^c 
+ ( X +7)3— 1 





AGm a 


( 1 + costf \ 




c 2 d a 





c 2 r ab \ 1 — sin $ 





2 G 2 m 2 / 1 + cos $ \ 2 G 2 m 2 sin$ (l + sin$ ) 



c 1 
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2 G 2 m a m\ ) sin $ cos i?o 



c A r ab d a l-sin^ 




(69) 



Similar to the case of 1-point-mass system, the 2PN light deflection in a 2-point-masses 
system with static approximation does not depend on PPN parameter j3, neither. One 
possible reason for this is that j3 usually associates with motion of bodies so that when 
the system could be treated as a static system due to T cross <C X, it would disappear in 
the deflection angle. However, when T croS s ~ X, we suppose the time evolutions of the 
gravitational fields, the motion of bodies and PPN parameter (3 would explicitly show up in 
the expression of deflection angle, and this more complicated issue will be the next move of 
our investigation. 

4- Quantitative examples 

If the light-ray just grazes the limb of the body a, the angle between light-ray 1 and 2 at 
the observer i? could be very close to and be neglected so that 




-(l + 7 ) 2(l+ 7 )--(7 + 8 7 )7r 




J cW a 



and it could be written in a more practical form 




(70) 



A# — 56qr + <5#stt, 



(71) 
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where 



4Gm a 15 M 2 m 2 a ^G 2 m a m h 

= L75 (S)(^:) arcsecond 

+3-50 ^V^V /.as 



M m y v d, 



a / 

2 



_ 2Gm a _ 31 N G 2 mg G 2 m a rnb , 2 . 

^ STT = 7 AT " 7(8 " " 107 rA^7 + ° (7 } 



2.1 x 10~ 5 )\M & )\d. 



— (^)te)fe) 2 (t)(^) - 

+C?(f). (73) 

If we ignore the mass of body 6, the results Eqs. ( 172]) and ( 1751) could return to the 1-point- 
mass case Eqs. (156"]) and ( 157j) . 

For experiments conducted in the Solar System, these results could be applied to some 
cases that only the Sun and the Jupiter need to be considered. If we suppose a light-ray 
passing the limb of the Sun, letting the body a be the Sun in Eqs. ( 172"]) and ( J75|) . then the 
2PN deflection angle of 5#gr contains 3.50 /xas caused by the Sun itself and —9.46 pico- 
arcsecond (pas) due to the coupling term of the Sun and the Jupiter. For the effect of the 
deviation from GR at the 2PN level in <5#stt, it consists of 81.1 pas from the Sun itself and 
—0.166 femto-arcsecond (fas) caused by the coupling effect of the Sun and the Jupiter, both 
of which are below the thresholds of future space missions. 

In another case that a light-ray just grazes the limb of the Jupiter, letting the body a be 
the Jupiter in Eqs. (172]) and ( 175]) . the deflection angle at 2PN level of 59gr consists of 0.302 
nas by the Jupiter and —0.0921 nas by the coupling of the Sun and the Jupiter. The 2PN 
contributions in 5#stt are made of 7.00 fas by the Jupiter and —1.61 fas by the coupling of 
the Sun and the Jupiter. 
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One of the major differences of this work from previous works is abandoning the hier- 
archized hypothesis, which leads to the 2-body effects in the light deflection at 2PN level. 
Although these effects in the 2PN light deflection is at least 1 order of magnitude lower than 
the thresholds of future ~nas experiments in the Solar System, their contributions in 5#gr 
and £#STT could reach —0.515 //as and —9.02 pas respectively for an imaginary observer 
in a binary system with two ~ 1 M Sun-like stars and separation by ~ 0.1 AU, both of 
which increase 3 orders of magnitude more than the system with mass ratio ~ 10~ 3 (such 
as Mj up itcr/M Q ) and the same separation. Hence, if this observer would be able to measure 
7 down to ~ 10~ 6 , the 2-body effect in 59 GR had to be considered, because the third term 
in £#gr and the first term in £#stt are comparable with the level of ~ 1 /ias. 

An even more extreme case would be replacing one Sun-like star in above hypothetical 
example with a neutron star with radius ~10 km, which could make the 2PN 2-body effects 
in 59gr and <5#stt raise to —35.9 mas and —0.627 //as. 

V. CONCLUSIONS AND DISCUSSION 

In this paper, we mainly focus on the scalar-tensor theory and obtain its 2PN metric. To 
investigate its next to the leading order contributions in the light propagation and related 
measurements in the near zone of a gravitational system, we reduce the metric to its 2PN 
form for light by supposing an A-point-masses system. In doing so, at 2PN level, we abandon 
the hierarchized hypothesis and do not assume two PPN parameters 7 and (3 to be unity, 
which makes our work differ from previous works. It is found that although there exist 7 
and (3 in the 2PN metric, only 7 appears in the 2PN equations of light. 

As one simple example of applications for future experiments, a gauge-invariant angle 
between the directions of two incoming photons for a differential measurement is investi- 
gated after the light trajectory is obtained in a static and spherically symmetric spacetime. 
It shows the deviation from GR 59 STT does not depend on f3 even at 2PN level in this 
circumstance, which is consistent with previous results. 

A more complicated application is light deflection in a 2-point-masses system. We con- 
sider a case that the light crossing time is much less than the time scale of its orbital motion, 
Across <S and thus treat it as a "static" system. The 2-body effect at 2PN level originating 
from relaxing the hierarchized hypothesis is calculated. Our analysis shows the 2PN 2-body 
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effect in the Solar System is one order of magnitude less than future ~ 1 nas experiments, 
while this effect could be comparable with 1PN component of 5#stt in a binary system with 
two Sun-like stars and separation by ~ 0.1 AU if an experiment would be able to measure 
7 down to ~ 10~ 6 . 

Our next move is to relax the static approximation for the 2-point-masses system with 
Across ~ X and investigate the time evolution pattern of the 2PN light deflection. We suppose 
the PPN parameter (5 would appear in the result. 
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Appendix A: Evolution equations of metric coefficients of second order post- 
Newtonian approximation 

In this appendix, we will give the metric coefficients of the scalar-tensor theory at the 
2PN order as follows. 
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(2) 

1. N and ( 

(2) 

The equations for N and ( are 



UN = -8nGa, 



(2) 

U( = -4(1 - j)irG<r, 
where □ is the D'Alembert operator in the Minkowski spacetime and 

= Up + 1 
coo + 2 

G 



0o(l + 7)' 
And we can see G = 1/0q when 7 = 1. 



2. Hij 



The equation for if yields 

UHij = SjirGcrSij, 
which gives = V5ij. Then, we obtain 

UV = SjirGa. 

Furthermore, we have V = 'yN by comparison with Eqs. (JSQ and dA6l) . 
3. Li 

The equation for is 

ULi = 8(7 + l)7rGa t . 



(4) 

4. L and £ 

(4) 

The equations for L and £ are 

/ 4/3_^_3(2)^ 
□L = -8ttGct \V - 2N + ^ — C ) - 8(7 - l)7rGa fcfe 
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where 



1 1 ( 2 ) 4(/3-l) ( 2 ) (2) 

~^N k N k - -N k V k - C , k N k - C , k ( , k , 



(4) 
□C 



= -AnGa 



8(7 - l)nGa, 



kk 



If If 4(1-/3) (2) (2) 

~*"2^ ' k ,k 2 ,fe (7 — l) 2 



1 + 



(2w + 3)(2w + 4)2- 



5. Q. 



y 



The equation for Qij reads as 



-8(1 + 7)7^ 



y 



-yy 



— NjN j + VNij — VV 

(2) 1(2) 1(2) 1(2) 1(2) 

- C JV y - - C - -CjN, + - c + - C jV, 

(2) (2) (2)(2) 2f2-y - 1)(2) (2) 
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4/3 -4 + 7 2 -7(2) 



+<Jij^ +8irGa 



+ ~f(N-2V) + 



7 
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1 1 (2) 4(5 - 1) (2) (2) 

+ ^N k V k + -V k V k - C , k Vk + -^7^ C, fc C,, !> . 



6. Qi 



The equation for Qi is 
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□Q. = +8(7 + l)7rG<Ti ^2V - JV - £ J + SttG^ 

+ W + ^iiV >t + + V t V t - 

— -^N )k L i)k + -N )ik L k — N )k L k)i + -V^fcLj jfc — -V^L* — V k L kj . 

1(2) 1(2) (2) 3(2) !(2) 

(2) (2) (3 7 _ 1)(2) (2) (2)(2) 

—Li >k C,fc — C ,ik^k h - _ j- C ,i C ,< + C C ,jt- 



27 



(6) 

7. Q and C 



(6) 

The equations for Q and ( are 



UQ = -8nGa{ + N 2 - 2NV - 2L - 



and 



4£- 

7 



(2) 

c 



+ 
+ 



1 - 



4(/3-l)( 7 2 + 8/3 - 2 7 - 7) 
( 7 -l)3 



-8nGa 



kk 



( 7 -l) 2 
4/3 -7 -3 ^ 
7-1 

7 2 , go _ 2 7 _ 7(2)- 

(2 7 -l)\/ + (2-7)iV- 7 " 7 - C 

7-1 



(2) Z 

c 



+167rG<T fc L fc - ^NN k N k + ^VN k V k 

+^N t N t + ^V t V t + N t V t + NN tt + VN tt 

—N^ k L k)t — L k N )kt — V kt L k — N ik L )k — -V k L tk 

1 

+N : hQki,i + N tk iQ k i — -NjQkkj — L k:l Li tk + L^ k L^ k 

(2) (2) (2) 

+2(, t V t -2( M L k - C, k L,k 

2( 7 2 -7 + 2/3-2) (y P) 4(/3 - 1) P) W W(2) 
+ ( 7 _ 1)2 C ,t C ,t + TTTTiy 7V C ,fc C ' fc C C ' fc ,fc 

26 WP) P) 4 (/3-l)( 7 2 + 8/3 - 2 7 - 7) ( 2 )( 2 ) ( 2 ) 
s s ,k s ,fc H 77 TvZ ^ s ,fc s ,fe 



+ 



(7 - 1) : 

(4) 



(7 " 1)^ 



S ,k S 



( 7 -l) s 



(6) 

□ C = -4ttGcx 



+ ( 7 - + L) - ?^-ii ( c - V) 

2(32/3 2 + i 7 - 64/3 - i + 32) ( 2 ) 2 + 8(/3 - 1) W 



—A^Ga kk 



(7 - I) 3 

16(5-1) ( 2 )' 
7 -l 3V-7V)- lP = J C 
7-1 

(2) (2) 



7-1 
- 87rG(7 - l)a k L k 



1 (2) 1 (2) (2) (2) 1 (2) 

+ 2^ C , k N k + -V C , k V k - 2 C M L k - C , fc L M + ^ C 

(2) (2) 1 (2) (2) 

+ C )fc Q fc i,i + C ,i k Qik - g C ,jQfcfc,j + (iv + v) c >tt 

4/3 + 2 7 - 7 2 - 5 ( 2 ) ( 2 ) 32(/3-l) 2 + 2i( 7 -l) ( 2 )( 2 ) ( 2 ) 
H ; — C ,t C ,t H ^ _ iy £ C ,fe C ,fc 



(7 - If 



(A13) 
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8(/3-l)(2) (4) 1(4) 1(4) 

JZ _ iy C ,fc C ,fc + 2 C ,k N ,k - 2 C ,k v ,k, 



(A14) 



where 



_ 1 (7 ~ I) 4 

L = W 2 . 

2 7 + 1 



(A15) 



In general relativity, i = 13 1. 



Appendix B: Application in a static, spherically symmetric spacetime 



1. Light-ray trajectory 



In this appendix, we will consider the gravitational field outside a static, spherically 
symmetric body and mainly focus on the parameters in 2PN level in several gravity theories. 
The trajectory of the light-ray can be obtained through integrating the following equation 

An a x)x 



x 



ft i \ Gm ° 

-( 1 + 7)-^r 



G 2 m 2 



2(1 + 7) + 



y n a x) 



[n n x x 



(Bl) 



by adopting the iterative method used in Ref. jl8(] . We assume the unperturbed light-ray 
as follows 



x N = x + c(t - t )k, 



(B2) 



where k l is a unit vector representing the light direction at t = — oo, t is an instant on the 
light path and x l Q is the position of photon at to- The instant to can be arbitrarily chosen, for 
example the emission or the observation instant. The photon's coordinates can be written 
as sum of perturbations on x l N \ 



x(t) = x N + 5x = x N + 5x 1PN + 5x 2 pn- 



And we use the following assumption for motion of the bodies 



(B3) 



r A (t) = x(t) -y a = x(t) - y a (t a ) 



(B4) 



where t a is the moment of the closest approach between the body a and the unperturbed 
light ray. y a (t a ) means the position of a-th body (y a ) is evaluated at the time t a . After 
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these, we obtain the results with the method as [18] 



-x(t) = k + -5x 1PN (x N ) + -5x 2PN (x N ) 

c c c 



(B5) 



x(t) = x N (t) + 



5x 1PN (x N ) - 5x 1PN (x ) 



+ 



5x 2PN (x N ) - 5x 2PN (x ) 



(B6) 



It is worthy of note that the 2PN terms in our solution actually have two sources, direct 
and indirect. The direct part comes from the 2PN order itself. The indirect part comes 
from the 1PN terms when the 1PN solution is iterated into itself in order to attain a 2PN 
accuracy, namely, we substitute xn + Sxipn into the trajectory of the light-ray in the 1PN 
approximation. And 
1 



-5x 1PN (x) 

c 

-5x 2PN (x) 

c 



gj.|fcx(r.xk) +fc 



c 2 r a 



k ■ r n 



(B7) 



1 G 2 m 2 G 2 m 2 

i ( i +7) £^( fc . r0r . + ^ (| . 1 (i +7 ) S 



>+7)(7 + 8 7 )l 



k ■ T n . 1 / 7T 

+ — — + arctan 



i ' d a \2 



a\' a 

k ■ r, 



2 

— + 

r„ 



1 



r a -k 



G 2 m 2 
+ ^^k 

err n 



-^(l + 7)(5 + 4 7 )--(l + 7) 2 - 
4 Tn r n 



k ■ r a 



(B8) 



and 



5xi PN (x) 
Sx 2PN (x) 



Gm a ( k x (r a x k) 
7)— —{ — r. + fcln(r a + k ■ r a ) 



cf 
G 2 m 2 



o( 1 + 7)-z-/r- a + 



r a -k-r a 

G 2 m 2 



(B9) 



c 4 r 2 



^k 



C L 



k ■ r n 



i(l + 7 )(7 + 8 7 )-3-arctan- 7 - 

o da da 



•(1+7)' 



k ■ r n 



k ■ r a y 



i(l + 7 )(7 + 87 )^ 



7r k ■ r a 

— h arctan 

2 d n 



(BIO) 



where d a = k x (r a x k) is an impact parameter to represent the closest distance between 
the unperturbed light ray and body a, d a = \d a \. For GR (7 = (3 = 1), Eqs. flB5j) and flB6j) 



will reduce to the results of Ref. 



id- 



2. Light deflection 



In some practical astronomical measurements, a differential measurement is more power- 

nn 

ful. This concept is employed by LATOR mission |1|, |2[ through a skinny triangle formed by 
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two spacecrafts and the International Space Station. Thus, we construct a gauge-invariant 



angle 9 between the directions of two incoming photons based on 

COS V = 



18 



4l|. It reads 



(Bll) 



where the spatial projection operator is 



h a /3 — 9 a/3 + U a Up, 



(B12) 



which projects the two incoming photons onto the hypersurface orthogonal to the observer's 
four- velocity u a = dx a /cdr and K° = dx±(t)/dt, and K 2 = dxf(t)/dt are the tangent 
vectors of the paths and x 2 (t) of the two incoming photons. Then, for a static observer, 
we obtain 

Gm a f fci • [fc 2 x (r x fc 2 )] . k 2 ■ [fci x (r„ x fci)] 



e = + (i + 7) 



c 2 r a sin i} 



k 2 ■ r a 



'(1 + 7) 5 



G 2 ml 
c 4 sin i?n 



k 2 • di a ( ^ ki ■ r a \ fci ■ d 2a ^ 



1 G 2 m 2 r 

~(l + 7) a 



c 4 r 2 



+-(7 + 8 7 )(l + 7) 

, ^2 • di a 



a la 

fci • r a k 2 ■ r a 

dia d 2a 

G 2 m 2 a (ki d 



2a 



fci ■ r a 
k 2 ■ r a 



2« 



d\ a 



7T 

— h arctan 
2 



c 4 sin t? 
fci • 



4, 



7T / fc 2 • r a 



— h arctan 
2 



2d 



C?la 



(B13) 



where 



arccos(fci ■ k 2 ) 



(B14) 



and we use that the position of the photon at the moment t of observation coincides with 
the position of the observer so that 



Xobs = Xqi + C(t - t l)fcl + SXi = X 02 + c(t - t 02 )k 2 + Sx 2 , 



(B15) 



where (toi^oi) denotes the moment and position of the light signal 1 of emission and 
(£02,^02) for the light signal 2 respectively. In Eq. ( 1B13I) . di a = fci x (r a x fci) and is an 
impact parameter to represent the closest distance between the unperturbed light ray 1 and 
a-th body. Similarly, d 2a = k 2 x (r a x k 2 ) and is an impact parameter to represent the 
closest distance between the unperturbed light ray 2 and a-th body. 
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Furthermore, we assume that one of the two light rays moves along the line connecting 
the body a and the observer, namely source 2. This means the impact parameter for source 
2 is zero so that d 2a = \k 2 x (r a x fc 2 )| = 0, k 2 ■ d la /d la = sin^ , k ± ■ r a /r a ~ cos^ , 
d\ a = r\ — (ki ■ r a ) 2 , ki ■ r a /d la w 7r/2 and di a /r a = 0. Then, we obtain Eq. (J51|) . 
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